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This is a short review of the quasinormal mode spectrum of Schwarzschild, Reissner-Nordstrom 
and Kerr black holes, summarizing results obtained in The summary includes previously 

unpublished calculations of i) the eigenvalues of spin-weighted spheroidal harmonics, and ii) quasi- 
normal frequencies of extremal Reissner-Nordstrom black holes. 
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Black holes are one of the most extreme predictions of Einstein's general relativity. They challenge common sense 
to such an extent that Einstein himself was reluctant to accept their reality. Despite the initial resistance, nowadays 
the existence of black holes is widely accepted by the physics community. Observationally, the astronomical evidence 
for black holes is compelling. Theoretically, a key step to accept the idea that black holes are not just a mathematical 
solution of the Einstein equations in vacuo was the proof of the (linear) stability of the Schwarzschild metric. This 
proof is based on black hole perturbation theory, as developed in two classic papers by Regge- Wheeler and Zerilli 0,1^. 

■ Their formalism gets rid of the angular dependence of the perturbation variables through a tensorial generalization of 
the spherical harmonics, thus reducing the solution of the perturbed Einstein equations to the study of a Schrodinger- 

• like wave equation with some potential V{r). Stability can then be established by methods that are familiar from 
^ \ quantum mechanics: perturbations due to an external (eg., electromagnetic or gravitational) field are considered as 

■ waves scattering off the corresponding potential |^ . It is worth stressing that the causally disconnected region inside 
I the event horizon plays no role whatsoever in this analysis. 

. The Regge- Wheeler /Zerilli formalism was later extended to encompass charged (Reissner-Nordstrom) or rotating 
I ' (Kerr) black hole solutions. For Reissner-Nordstrom (RN) black holes we have two coupled Schrodinger equations, 
, describing the mutual interaction of gravitational and electromagnetic perturbations. For Kerr black holes the problem 
' is highly non-trivial, since the background geometry is not spherical and tensor spherical harmonics cannot be used 
^ \ any more. Nonetheless, Teukolsky was able to separate the angular and radial dependence of the perturbations using 
O^" a certain class of special functions (the spin-weighted spheroidal harmonics; see section IV C|l - a truly remarkable 
\ achievement • A complete discussion of black hole perturbation theory can be found in the monumental monograph 
bJ[)i by Chandrasekhar I3| . The extension to charged and rotating (Kerr-Newman) black holes turned out to be a formidable 
task: the numerous attempts to separate the corresponding perturbation equations have all failed. 

Once we know that a black hole solution is indeed stable, the next step is to understand "how much stable" it is. 
In other words, on which timescale does a black hole radiate away its "hair" after formation? The answer to this 
question is provided by the black hole's quasinormal mode spectrum. After an initial phase depending on the details 
of the collapse, the black hole starts vibrating into "quasinormal" (exponentially decaying) oscillation modes whose 
frequencies and decay times depend only on the intrinsic features of the black hole itself, being insensitive to the 
details of the collapse. 

We have seen that for Schwarzschild, RN and Kerr black holes the perturbation problem can be reduced to the study 
of certain characteristic wave equations. The background spacetime curvature determines the shape of the potential 
in each wave equation. Quasinormal modes are defined as solutions that satisfy "natural" boundary conditions for 
classical waves: they are purely outgoing at spatial infinity and purely ingoing at the black hole horizon. In the Fourier 
domain, this is an eigenvalue problem for the wave's (complex) frequency lo = w^j -I- iw/ . The two boundary conditions 
are satisfied by a discrete set of quasinormal frequencies, uj = {wn}, where the integer n = 0, 1, . . . labels the modes, 
which are normally sorted by increasing values of Quasinormal frequencies are complex {quasi-novvcial) because 
outgoing waves at infinity radiate energy (the system is not conservative). The damping time of each oscillation is 
of course given by t = l/|a;/|: modes having large values of jw/j damp faster. The quasinormal mode spectrum, 
unlike ordinary normal mode expansions, is not complete: in general, we cannot express a given perturbation as a 
superposition of quasinormal modes (see section 4.4 of reference Q for a discussion). Completeness would be desirable 
from a mathematical point of view, but in any case quasinormal mode expansions are useful (and routinely used) 
to describe the final stages of formation of a black hole after collapse. For example, it has been shown that the 
numerical waveforms obtained from general relativistic simulations of distorted black holes are very well fitted using 
the fundamental quasinormal mode, with the possible addition of a few overtones p^ . 

The rest of this short review is devoted to a detailed presentation of the quasinormal mode spectrum of 
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Schwarzschild, RN and Kerr black holes. The quasinormal spectrum has a complicated structure: as it happens, 
black hole perturbation theory leads us once again "into a realm of the rococo: splendorous, joyful, and immensely 
ornate" 01 ■ Such an ornate spectrum can at times be hard to compute, even by the most sophisticated numerical 
methods. 

In particular, for a long time the high-damping limit of the spectrum has been a no-go zone. From a physical 
point of view, it is clear that slowly-damped quasinormal modes should dominate the response of a black hole to 
any type of perturbation: this obvious consideration initially reduced the impetus to explore the large- Iw/I (high 
damping) regime. From a technical perspective, most numerical methods to impose the quasinormal mode boundary 
conditions are based on discriminating outgoing waves from ingoing waves. The exponentially decaying solution is 
quickly swamped into numerical noise, and this effect is larger for larger values of the damping. As a consequence, 
all numerical methods are doomed to failure at some critical value of \oJi\. 

More recently, the high-damping region of the quasinormal mode spectrum gained more attention because of a 
conjectured connection with quantum gravity. York first attempted to relate the (purely classical) quasinormal 
mode spectrum with quantum properties of black holes, and more specifically Hawking radiation. In 1998 Hod 
noticed that, as jw/j — > oo, the numerically computed oscillation frequencies of Schwarzschild black holes tend to a 
constant asymptotic value wjj = Th In 3, where Th is the Hawking temperature of the black hole. In the same limit, 
the (constant) spacing between modes is given by 2ttTh- Assuming that the minimum energy emitted or absorbed by 
a black hole is AM = fiw^ and using the Bekenstein-Hawking area-entropy relationship. Hod found an equally spaced 
black hole area spectrum (in agreement with some heuristic arguments by Bekenstein 0|). By virtue of the "In 3" 
factor, the resulting spectrum is also consistent with a statistical-mechanical interpretation of black hole entropy. The 
idea was revived more recently by Dreyer [14| , who used a similar argument to fix the Barbero-Immirzi parameter [T^ 
- an undetermined constant of proportionality entering black hole entropy calculations in Loop Quantum Gravity. 

The results we show here do not seem to support an extension of Hod's original proposal to charged and rotating 
black holes. In addition, the original Loop Quantum Gravity calculation of the black hole area spectrum was recently 
re-examined and shown to be in error by Domagala and Lewandowski These developments could be used to 

argue against the relevance of quasinormal modes for Loop Quantum Gravity, or any other theory of quantum gravity 
for that matter. Some hints to clarify - or kill once and for all - the connection between quasinormal modes and 
area quantization might come from analogue black holes |16| , but in that case the thermodynamical interpretation of 
horizon area as entropy is not a trivial issue. 

Whether or not quasinormal frequencies have something to say about quantum gravity, the recent conjectures 
stimulated a lot of work to clarify the structure of the high-damping region of the quasinormal mode spectrum. The 
following sections summarize some of these recent developments. Our aim is not to give a complete overview of 
recent work in this field: that would be an ambitious goal. For example, we do not discuss recent progress on the 
asymptotic quasinormal spectrum of black holes in non-asymptotically fiat spacetimes or higher dimensions (with a 
few exceptions where appropriate). We only present a personal (and therefore biased) account of some open problems, 
concentrating on the Kerr-Newman family of black hole solutions. Before turning to a description of the quasinormal 
mode spectrum, it is useful to mention in passing different numerical methods that have been used in the past for 
quasinormal mode calculations. 



The summary we present in this section is only meant to highlight pros and cons of some popular numerical methods 
to compute quasinormal modes, and ultimately to justify our use of the continued fraction technique in the present 
context. The list is by no means exhaustive. For a more complete discussion of the merits and drawbacks of each 



1) Time-domain evolutions of the Rcgge- Wheeler and Zerilli equations were first performed by Vishveshwara [Tg 
and Press |il9i]: this is actually how quasinormal modes were first discovered. However, finding highly damped 
modes in this way requires long and stable evolutions. This method is not used any more for quasinormal 
mode calculations, with some notable exceptions in cases where a Fourier-domain treatment is cumbersome or 
impossible. Examples include i) the simulation of complicated distributions of matter in a fixed (black hole 
or stellar) background and ii) time evolutions of the Kerr perturbation equations (see [IJ and references 
therein). The latter are particularly important for the simulation of extreme mass ratio inspirals, one of the 
target sources for the space-based gravitational wave interferometer LISA. No generalization of the Carter 
constant is known for off-equatorial, non-circular particle orbits in the Kerr metric, and time evolutions provide 
a valid alternative to calculations of radiation reaction effects using self-force methods. 



II. NUMERICAL METHODS 
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2) Direct integrations of the wave equation in the frequency domain were pioneered by Chandrasekhar and Detweiler 
|22|. The simplest version of this approach is numerically unstable: the exponentially vanishing solution is very 
hard to extract from numerical noise in the exponentially growing solution, a problem that plagues especially 
highly damped modes. Several improvements have been attempted. A notable one was introduced by NoUert 
and Schmidt, who used a sophisticated Laplace transform definition of the quasinormal frequencies [23]. In 
this framework, the state of the art is probably the complex integration technique developed by Andersson and 
coworkers [23 |. In any case, direct integrations are quite time consuming. Phase-integral based methods (see 
below) sometimes involve numerical integrations; they can be competitive with continued fraction methods at 
moderate dampings, but ultimately fail at very high values of 

3) Inverse potential methods have been introduced by Mashoon and collaborators |25j |. The idea is to approximate 
the exact potential by some simpler potential (eg. a potential of the Poschl- Teller form) whose spectrum can 
be found analytically. These techniques can be a valuable tool in certain situations 26J, but in general results 
obtained from approximating potentials should not be expected to be accurate. For example, NoUert showed 
that approximating the Regge- Wheeler potential by a series of square wells of decreasing width does not yield 
the expected quasinormal frequencies as the width of the square wells tends to zero. Analytic approximations 
can still be useful in certain limits, such as the eikonal approximation (large-Z limit) or the high damping limit 



4) WKB methods were originally suggested by Schutz and Will |3C| and are based on elementary quantum mechani- 
cal arguments. The analogy between the wave equation describing black hole perturbations and the Schrodinger 
equation can be exploited to use standard WKB expansions. At lowest order, quasinormal frequencies are given 
by the Bohr-Sommcrfcld rule: 



where rA and rn are the two roots (turning points) of — V{r) — 0. The WKB expansion can be pushed to 
higher orders [S^, |32| , allowing a simple and reasonably accurate determination of the quasinormal frequencies 
and a trasparent physical interpretation of the results. The method can be generalized to deal with RN [s^ 
and even Kerr black holes js^]. In the latter case, unfortunately, the WKB method is not very simple to apply 
(the Kerr potential in the standard Teukolsky equation is complex in the first place, but it can be made real 
transforming the equation to a different form) and does not yield very accurate results. The WKB approximation 
gets better for large values of I, allowing an analytic determination of quasinormal frequencies in the large-l 
limit, but it was found to fail badly at high values of the damping js^. 

5) Phase-integral methods are a substantial improvement over standard WKB techniques |36ll37ll3^. Their major 
drawback is related to the fact that they require integrations in the complex plane. The integration contour 
must be found case by case, and the problem can be quite tricky, depending on the singularity structure of the 
potential. A variant of the phase-integral method has been used to compute analytically the Schwarzschild and 
RN quasinormal frequency in the limit jtj/l — > oo [s^. A similar analytical treatment for the Kerr case is still 
missing. 

6) Continued fractions were first used by Leaver |4Clj | exploiting once again an analogy with quantum mechanics. 
Leaver's approach is based on a classical 1934 paper by Jaffe on the electronic spectra of the 11+ ion [4ll| . 
The key steps are: a) recognize that the wave equations can be seen as special cases of generalized spheroidal 
wave equations; b) write appropriate series expansions for the solutions; c) replace the series expansions in 
the differential equations and derive recursion relations for the expansion coefhcients; d) use the recursion 
relations to study the convergence properties of these series and impose quasinormal mode boundary conditions. 
The determination of quasinormal frequencies boils down to a numerical solution of continued fraction relations 
involving the mode frequency and the black hole parameters. The evaluation of continued fractions only involves 
elementary algebraic operations and the convergence of the method is excellent, even at high damping. In the 
Schwarzschild case, in particular, the method can be tweaked to allow the determination of modes of order up to 
~ 100, OOP [4^ . Gauss- Jordan elimination techniques can be used to encompass the RN 0| and even extreme 
RN lH cases. 

Other approaches may be preferred in specific situations, but Leaver's continued fraction technique is the best 
workhorse method to compute highly damped modes in general situations. It is very reliable, largely independent 
of the form of the potential (that is only used to find the recursion coefficients once and for all) and fast, bypassing 
numerical integrations of the wave equations. 




(1) 
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All the results we present in the following have been obtained using different variants of Leaver's method. We 
adopt geometrical units and consistently use Leaver's conventions piOi] : in particular, we set 2M = 1. This means that 
extremal Kerr and RN black holes correspond, respectively, to a = 1/2 and (5 = 1/2. 



III. SCHWARZSCHILD BLACK HOLES 



A. Computational method 

Let us begin from the simplest case of spherically symmetric, uncharged, non-rotating Schwarzschild black holes. A 
detailed derivation of the perturbation equations can be found in j^J (for more details on the computational procedure 
see ^3,E3L-I^ brief: the angular dependence of the metric perturbations can be separated using tensorial spherical 
harmonics |5| . Depending on their behavior under parity, the perturbation variables are classified as polar (even) or 
axial (odd). We can get rid of the time dependence using a standard Fourier decomposition for each perturbation 
variable ^p^. 

1 /■+°° 

= e-''^'^K^,r)duj. (2) 

The resulting differential equations can be manipulated to yield two wave equations, one for the polar perturbations 
(that we shall denote by a superscript plus) and one for the axial perturbations (superscript minus): 

£^^+.^'^Z^^V^Z\ (3) 

where we introduced a tortoise coordinate defined in the usual way by the relation 

— (4) 

and A = r(r — 1). 

Polar perturbations are related to axial perturbations by a differential transformation discovered by Chandrasekhar 
0. As anticipated in the introduction, quasinormal modes are solutions of equation (|3Jl that are purely outgoing at 
spatial infinity (r —^ oo) and purely ingoing at the black hole horizon (r — > 1). These boundary conditions arc only 
satisfied by a discrete set of complex frequencies (the quasinormal frequencies). For the Schwarzschild solution the 
two potentials are quite different, yet the quasinormal modes for polar and axial perturbations are the same j^. 
This can be seen as a consequence of the two potentials being related by a supersymmetry transformation [45l | . Since 
polar and axial perturbations are isospectral and V~ has an analytic expression which is simpler to handle, we can 
concentrate on the axial equation. Written in terms of the "standard" radial variable r, this equation reads: 



r(r - 1)—^ + — 



l{l + l)- 



1 



= , (5) 



where s is the spin weight of the perturbing field {s — 0, — 1, — 2 for scalar, electromagnetic and gravitational 
perturbations, respectively) and I is the angular index of the perturbation. Notice that perturbations of a Schwarzschild 
background are independent of the azimuthal quantum number m, because the background spacetime has spherical 
symmetry; the same is true in the RN case, but not for Kerr black holes. Equation Q can be solved using a series 
expansion of the form: 

4>i = ir- i)-i-r2i^e'"('-i) J2 «J ( ) • (6) 

j=o V ^ / 

where the prefactor is chosen to incorporate the quasinormal boundary conditions at the two boundaries. Substituting 
the series expansion (j^ in (jSJ we get a three term recursion relation for the expansion coefficients Qj : 

aofli + PqQq = 0, (7) 
ajaj+i + (3jaj + ^jUj-i =0, j = 1, 2, . . . 
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FIG. 1: Quasinormal frequencies for gravitational perturbations with Z = 2 (blue circles) and Z = 3 (red diamonds). Compare 
eg. Figure 1 in 42]. In both cases we mark by an arrow the algebraically special mode, that is given analytically by Equation 
I13I I: a more extensive discussion of this mode is given in section IVTI Notice that as the imaginary part of the frequency tends 
to infinity the real part tends to a finite, I -independent limit. 



where aj, (5j and 7j are simple functions of the frequency I and s. Their explicit form can be found in |4Q|. The 
quasinormal mode boundary conditions are satisfied when the following continued-fraction condition on the recursion 
coefficients holds: 

. ao7i "172 

The n-th quasinormal frequency is (numerically) the most stable root of the n-th inversion of the continued-fraction 
relation i.e., it is the root of 

a an-lln OLn-2ln-\ ^071 an7n-|-l an-|-l7n+2 I ^ ^ \ /n\ 

-r. ■■■^ = -n n ■•■ (n = I,2,...). (9) 

Pn-1— Pn-1— PO Pn+1— Pn+2 — 

The infinite continued fraction appearing in equation Q can be summed "bottom to top" starting from some large 
truncation index N. NoUert has shown that the convergence of the procedure improves if the sum is started using 
a wise choice for the value of the "rest" of the continued fraction, Rn, defined by the relation 

R. = """^ ^ . (10) 

PN+l — aN+l-tlN+l 

Assuming that the rest can be expanded in a series of the form 

oc 

i?Ar = ^CfciV-'=/2^ (11) 

k=0 



it turns out that the first few coefficients in the series are Cq = —1, Ci — ±\/—2iuj, C2 = (3/4 + 2iw) and C3 — 
+ l)/2 + 2uj^ + 3i(jj/2 + 3/32] /Ci (the latter coefficient contains a typo in ^42], but it is numerically irrelevant 
anyway) . 



B. Results 



We computed quasinormal frequencies using Leaver's technique as improved by NoUert. The slowly damped modes 
of the resulting spectrum for I = 2 and I = 3 are shown in Table and Figure ^ Reintroducing physical units, the 
fundamental oscillation frequency / = Wi{/(27r) and the damping time r = of an astrophysical black hole scale 

with mass according to the relation 

1^207 (i2ife) kHz, 0,5537 (^)„„^ (12) 



M J ' V 10 M( 
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TABLE I: Representative Schwarzschild quasinormal frequencies for / — 2 and / = 3 (from (4^). 





/ = 


= 2 




/ = 3 




n 












1 


(0 


747343, 


-0.177925) 


(1.198887, 


-0.185406) 


2 


(0 


693422, 


-0.547830) 


(1.165288, 


-0.562596) 


3 


(0 


602107, 


-0.956554) 


(1.103370, 


-0.958186) 


4 


(0 


503010, 


-1.410296) 


(1.023924, 


-1.380674) 


5 


(0 


415029, 


-1.893690) 


(0.940348, 


-1.831299) 


6 


(0 


338599, 


-2.391216) 


(0.862773, 


-2.304303) 


7 


(0 


266505, 


-2.895822) 


(0.795319, 


-2.791824) 


8 


(0 


185617, 


-3.407676) 


(0.737985, 


-3.287689) 


9 


(0 


000000, 


-3.998000) 


(0.689237, 


-3.788066) 


10 


(0 


126527, 


-4.605289) 


(0.647366, 


-4.290798) 


11 


(0 


153107, 


-5.121653) 


(0.610922, 


-4.794709) 


12 


(0 


165196, 


-5.630885) 


(0.578768, 


-5.299159) 




/n 
(U 


i looua, 


-y.DDUo r y ) 


(U.4U4io r , 


n QQQ101^ 

-y.oooizi ) 


30 


(0 


165814, 


-14.677118) 


(0.257431, 


-14.363580) 


40 


(0 


156368, 


-19.684873) 


(0.075298, 


-19.415545) 


41 


(0 


154912, 


-20.188298) 


(-0.000259 


,-20.015653) 


42 


(0 


156392, 


-20.685630) 


(0.017662, 


-20.566075) 


50 


(0 


151216, 


-24.693716) 


(0.134153, 


-24.119329) 


60 


(0 


148484, 


-29.696417) 


(0.163614, 


-29.135345) 



An "algebraically special" mode, whose frequency is (almost) pure imaginary, separates the lower quasinormal mode 
branch from the upper branch 46] . This algebraically special mode has a very peculiar nature, as we will see in 
section IVll It is located at 

- ,. (/-l)/(/ + !)(? + 2) 

f^i = ±i Q , (13) 

and it can be taken as roughly marking the onset of the asymptotic high-damping regime. The algebraically special 
mode quickly moves downwards in the complex plane as I increases: from Table n we see that it corresponds to an 
overtone index n — 9 when 1 = 2, and to an overtone index n = 41 when Z = 3. This means that for high values of I 
the asymptotic high-damping regime sets in later, becoming harder to probe using numerical methods. 

NoUert was able to compute for the first time highly damped quasinormal frequencies corresponding to gravitational 
perturbations 42]. His main result was that the real parts of the quasinormal frequencies are well fitted, for large n, 
by a relation of the form 

UJR = LOoo + ■ (14) 

The leading-order fitting coefiicient ujao = 0.0874247 is independent of /. Corrections of order ~ n~^^'^, however, are 
Z-dependent (we will see in a moment that they also depend on the spin s of the perturbing field; that's why we 
denoted them by As,;). For gravitational perturbations (s = —2) NoUert found that X-2,2 = 0.4850, A_2,3 = 1.067, 
A_2,6 = 3.97. Numerical data for scalar perturbations (s = 0) are also well fitted by formula (|14() . and the asymptotic 
frequency for scalar modes is exactly the same. Only the leading-order correction coefficients Aq,; are different: 
Ao,o = 0.0970, Aoa = 0.679, Ao,2 = 1-85 

These numerical results are perfectly consistent with analytical calculations. Motl analyzed the continued 
fraction condition ((SJ to find that highly damped quasinormal frequencies satisfy the relation 

w - In 3 + (2n + l)mTH + OirT^I'^) . (15) 

(remember that in our units 2M = 1, so that the Hawking temperature of a Schwarzschild black hole Th = l/47r). This 
conclusion was later confirmed by complex-integration techniques and phase-integral methods |39j ]. Neitzke p9| 
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first suggested that leading order corrections to the asymptotic frequency should be proportional to [{s"^ — !) — 31(1 + 1)]; 
this is consistent with the numerical values of the Xg/s. The proportionality constant is dependent on s: A^,; = 
ks[{s'^ — 1) — 31(1 + 1)]. Maassen van den Brink ji^ derived k^2 analytically, finding 



V2[r(i/4)]^ 



-0.0323356. 



(16) 



4327r5/2 

in perfect agreement with the numerical data, that yield /c_2 — —0.0323, ko — —0.0970 = 3fc_2 "ll- 

Finally, a systematic perturbative approach to determine lower-order corrections to formula H14(l has been developed 
in [4^ . Their final result, valid for scalar and gravitational perturbations, is again perfectly consistent with numerical 
data. The perturbative tec hniq ue used in (47j has also been extended to higher-dimensional Schwarzschild-Tanghcrlini 
black holes with the result 4«| 



OJ 



-{D~Z)/{D-2) 



(17) 



Th ln(l -I- 2 cos TTj) -I- {2n -f 1)^7// + koLOi 

where the coefficient ku can be determined analytically for given values of D and j '49*1 . 

Leaver's method can easily be generalized to these higher-dimensional black holes i48ih the recursion relation has 
four terms for vector-gravitational and tensor-gravitational perturbations (see also [50j for similar calculations using 
WKB methods). The most impressive example of an application of Leaver's technique to date is probably reference 
|5lj | . There the continued fraction method has been used for a detailed calculation of scalar-gravitational quasinormal 
frequencies in the case of 5-dimensional Schwarzschild black holes using an eight-term recursion relation. 

In conclusion, all numerical and analytical results for non-rotating black holes are in perfect agreement. As |a;/| 00 
scalar and gravitational oscillation frequencies of non-rotating Schwarzschild- Tangherlini black holes in all dimensions 
tend to the constant value ojr = Th In 3 (this was first argued in ||23] and explicitly shown in [^). For electromagnetic 
perturbations the situation is less clear: analytical and numerical results in D = 4 suggest that the asymptotic 
limit should be uir — 0. For electromagnetic perturbations in = 5, formula H17|l doesn't even make sense: this 
probably means that there are no asymptotic quasinormal frequencies at all (a possibility first suggested in "2^). A 
similar behavior occurs for the 3-dimensional "draining bathtub" metric describing a rotating acoustic black hole 16], 
suggesting that this feature may somehow be related to odd-dimensional spacetimes. 



IV. REISSNER-NORDSTROM BLACK HOLES 



A. Computational method 



The computational procedure described in the previous section can be generalized to RN black holes. Here we only 
sketch the key steps, referring the reader to 43] for more details. As usual we define a tortoise coordinate by the 



relation |0J, but now A 

be written as 



r + Q (recall that in our units < Q < 1/2). Explicitly, the tortoise coordinate can 



■ ln(r — r+) — 



■ ln(r — r_) 



(18) 



where r± — {1 ± \/l — 'iQ^)/2 is the location of the inner (Cauchy) and outer (event) horizons of the RN met- 
ric. After separation of the angular dependence and Fourier decomposition, electromagnetic and axial-gravitational 
perturbations of the RN metric are described by two wave equations. 



where 



and 



V- 



A 



dr'i_ 



l{l + l)r - qj + 



Vrzr 



4Q2 



i,j = l,2 (i^j) 



qi = 



3+ v/9 + 16Q2(^_i)(^ + 2) /2 



92 



3- v/9 + 16Q2(^_i)(^ + 2) /2 



(19) 



(20) 



(21) 



As in the Schwarzschild case, the polar perturbation variables (Z^) can be obtained from the axial variables (Z- ) by 
a Chandrasekhar transformation [g, so we do not consider them in the following. In the limit Q — > the potentials 
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V-[' and describe, respectively, purely electromagnetic [spin weight s = — 1 in Equation ISJ] and pure axial- 
gravitational [spin weight s = —2 in Equation (jSJ] perturbations of a Schwarzschild black hole; but for any non-zero 
value of the charge, electromagnetic and gravitational perturbations are tangled. The perturbation equations are 
solved by a series expansion of the form 

Zr = ^ ^_u~^^rl/ir^-r-)J2a,u\ (22) 

^ 3=0 

where u — {r — r^^) / [r — r ^) . The coefficients aj of the expansion are now determined by a four-term recursion 
relation whose coefficients depend on the angular index Z, on the frequency oj and on the charge Q. The problem can 
be reduced to a three-term recursion relation of the form (jSJ using a Gaussian elimination step (see |43(| for details). 
Then we can use again a continued- fraction condition to determine quasinormal frequencies, and NoUert's procedure 
can be applied to accelerate convergence for large values of the damping. The first few coefficients in the series (jllf) 
are now Co = -1, Ci = ±^-2\uj{2r+ - 1), C2 = (3/4 4- 2iwr+). 

Due to convergence reasons, continued fraction calculations become more and more computationally intensive for 
large values of the charge f43|. In the limit of maximally charged black holes {Q = 1/2) the wave equations have a 
different singularity structure, and deserve a special treatment j4^. We discuss the peculiarities of the extremal RN 
case in section HV CI 



B. Results 



The first few overtones of a RN black hole were studied numerically by Andersson and Onozawa |33, |53| using 
different numerical methods 54]. These studies revealed a very peculiar behavior. 
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FIG. 2: Trajectory described in tlie complex-frequency plane by the fundamental RN quasinormal mode as the charge is 
increased. The solid line corresponds to Z = 2 and Z2 (perturbations that reduce to the axial-gravitational Schwarzschild case 
as Q ^ 0); the dashed line, to Z = 1 and (purely electromagnetic Schwarzschild perturbations in the limit Q ^ 0). The 
modes coalesce in the extremal limit Q 1/2. 



The solid line in Figure|21is the trajectory described in the complex-frequency plane by the fundamental quasinormal 
mode with I — 2 corresponding to (perturbations that reduce to the axial-gravitational Schwarzschild case). The 
dashed line is the same trajectory for the fundamental quasinormal mode with I = 1 corresponding to Z^ (which 
limits us to purely electromagnetic perturbations). The Schwarzschild limit corresponds to the bottom left of each 
curve, and the trajectories are described counterclockwise as Q increases. The real part of the frequency grows 
monotonically with Q, and the imaginary part shows a relative minimum. This behavior can easily be understood 
using WKB arguments: restoring the mass in the equations, the lowest order WKB approximation yields 
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FIG. 3: The top two panels show the behavior of the n = 5 and n — 10 quasinormal frequencies in the complex lu plane. The 
n = 10 mode "spirals in" towards its value in the extremal charge limit; the number of spirals described by each mode increases 
roughly as the mode order n. The panels in the second row show how the n = 10 spiral "unwinds" as the angular index I is 
increased (in other words, the asymptotic behavior sets in later for larger I's). In all cases, we have marked by an arrow the 
frequency corresponding to the Schwarzschild limit (Q = 0). 



where Tq ~ (3M + •\/9M^~^^"8Q^) /2 is the location of the potential peak. It is worth noting that modes of with 
angular index I coalesce with modes of with index (l — 1) in the extremal limit. This is a general feature we are 
going to discuss in section ITV CI 

The algebraically special frequency of Reissner-Nordstrom black holes can be found using the same reasoning as in 
the Schwarzschild case j42j. It reads 



^ ^. {l^mi+l)il + 2) 



l±v/l + 4Q2(/ -!)(/ + 2) 



(24) 



where the superscript (i) refers to the two different potentials in Equation H19() . For extremal black holes {Q — 1/2) 
the first algebraically special mode ioi l — 2 is located at ft = — 3i. 

In the rest of this section we explore the high-damping regime. We fix our attention on the wave equation for Z2 
(perturbations reducing to gravitational perturbations of Schwarzschild as Q — s- 0); results for Z^ are similar. As can 
be seen in Figures |3| El and |S1 the behavior of high order modes is quite surprising. At first the trajectories described 
by the modes in the complex-cj plane show closed loops, as in the top left panel of Figure|21 Then they get a spiral-like 
shape, moving out of their Schwarzschild value and looping in towards some limiting frequency as Q tends to the 
extremal value. This kind of behavior is shown in the top right panel of Figure O We observe that such a spiralling 
behavior sets in for larger values of the modes' imaginary part (i.e., larger values of n) as the angular index / increases. 
In other words, increasing l for a given value of the mode index n has the effect of unwinding the spirals. This can be 
seen in the two bottom panels of Figure 13 However, for each / the spiralling behavior is eventually observed when n 
is large enough. 
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FIG. 4: Real part of the RN quasinormal frequencies as a function of charge for n = 5, 10, 30, 60, 5000, 10000, 100000. As 
the mode order increases the computation becomes more and more time consuming, the oscillations become faster, and a good 
numerical sampling is rather difficult to achieve; therefore in the last plot we use different symbols (small squares, circles and 
triangles) to display the actually computed points. For n = 5000, 10000, 100000 we also compare to the prediction of the 
analytic formula 12511 derived by Motl and Neitzke 28] . The oscillatory behavior is reproduced extremely well by their formula, 
but the disagreement increases for small charge: formula 12511 does not yield Th In 3 in the Schwarzschild limit. 



A perhaps clearer picture of the modes' behavior can be obtained looking separately at the real and imaginary parts 
of the mode frequencies as function of charge. Let us focus first on the real parts. The corresponding numerical results 
are shown in Figure 0| It is quite apparent that, as the mode order grows, the oscillating behavior as a function of 
charge start earlier and earlier. For large n the oscillations become faster, the convergence of the continued fraction 
method slower, and the required computing time gets longer. Therefore, when the imaginary part increases it becomes 
more difficult to follow the roots numerically as we approach the extremal value Q — 1/2. That's why data for very 
large values of n do not cover the whole range of allowed values for Q. Despite these difficulties, we have many reasons 
to trust our numerics. We have carefully checked our results, using first double and then quadrupole precision in 
our Fortran codes (indeed, as n increases, we can obtain results for large values of the charge only using quadrupole 
precision). Our frequencies accurately reproduce NoUert's results in the Schwarzschild limit, so the numerics can 
be trusted for small values of Q. More importantly, at large values of the charge numerical data are in excellent 
agreement with analytical predictions "2^. Motl and Neitzke found the following result for the asymptotic frequencies 
of scalar and electromagnetic-gravitational perturbations of a RN black hole |53 : 

e'^" + 2 + ie-'^"^ = 0. (25) 

This implicit formula for the asymptotic quasinormal modes was later rederived by an indepent method in !39| . The 
complex solutions of Equation H25|l exactly overlap with the oscillations we observe at large damping for large enough 
values of Q (see Figures 0] and (S)) . This agreement gives support to the asymptotic formula H25|l and confirms that 
the numerical calculation is still accurate for large values of Q 29] . 

Similar considerations apply to the imaginary parts. Some plots illustrating the general trend are shown in Figure 
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FIG. 5: Imaginary part of the RN quasinormal frequencies as a function of charge for n = 10, 30, 60, 5000. For n = 5000 we 
also display the actually computed points, and compare to the prediction of the analytic formula l|25|l . As for the real part, 
the oscillations are reproduced extremely well, but the disagreement with our numerical data increases for small charge. 



Once again the numerics show excellent agreement with the asymptotic formula H25|l as n increases. The analytic 
formula deviates from numerical results only for small values of the charge. 

In conclusion, numerical results lend strong support to the analytical expression (|25|l . This formula presents us 
with some striking puzzles. As emphasized in [l^llsl, the predicted asymptotic RN quasinormal frequencies do not 
reduce to the expected Schwarzschild limit. Equation (|15|l . One finds instead 



Tff ln5 



as 



(26) 



a prediction that is hard to reconcile with Hod's interpretation of the asymptotic quasinormal spectrum. Surpris- 
ingly, in the extremal limit Q — > 1/2 the real part of the frequency predicted by Equation H25|l coincides with the 
Schwarzschild value: 



TfflnS as Q^l/2. 



(27) 



It is not clear to which extent this result is relevant: in the extremal limit the two horizons coalesce, the topology of 
the singularities in the differential equation changes, and the problem may require a separate analysis. Indeed, there 
are some arguments that the asymptotic oscillation frequency for extremal RN black holes is not given by Equation 

Another bizarre prediction is that - according to Equation H25|l - asymptotic quasinormal frequencies of charged 
black holes depend not only on the black hole's Hawking temperature, but also on the Hawking temperature of the 
(causally disconnected) inner horizon. The dependence on the causally disconnected region is not surprising when 
one looks at the calculation by Motl and Neitzke: their result depends only on the behavior of the potential close to 
the curvature singularity. Still, the region inside the event horizon plays no role whatsoever in the classical analysis 
of black hole perturbations. From this perspective, the appearance of /3/ in Equation (|25|l is somehow disturbing. 
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FIG. 6: Quasinormal frequencies for gravitational-type perturbations of an extremal RN black hole with I = 2 (blue circles) 
and / = 3 (red diamonds). The spectrum for electromagnetic-type perturbations with I — 1 and / = 2 would be exactly the 
same. We mark by an arrow the location of the algebraically special mode. We find a quasinormal mode at tj = (0, —3.047876) 
[filled blue circle], while Chandrasekhar's formula (|24ll predicts a mode at a; = (0,-3) [empty blue circle]. Unfortunately the 
numerical method becomes unstable for values of \uji\ > 10, and cannot be used to verify the prediction (1271 of the Motl-Neitzke 
formula (I25II . 



An interesting classification of the solutions of Equation (|25|) can be found in Figure 3 and section 4.4 of [3^. Besides 
"spiralling" solutions the equation also admits periodic solutions when k = -\/l^^Q^/M^ is a rational number, and 
even pure imaginary solutions that may not be quasinormal modes at all! 

To conclude this section: analytical and numerical calculations indicate that (|25|1 is correct. Still, we are far from 
a full understanding of the physical meaning of this equation and its relation with the asymptotic Schwarzschild 
quasinormal spectrum (|15|l . 

C. Extremal Reissner-Nordstrom black holes 

Leaver's method fails in the extremal limit. In this limit the two horizons coalesce (rhor — = r_ = 1/2). The 
tortoise coordinate becomes, in our units, 

n =r + ln(2r-l)-^-^, (28) 

and the radial wave equation has irregular singular points at the horizon and at infinity. The series expansion l|22(l 
is not valid any more. The trick used by Onozawa et al. |44| is to change the expansion variable: they expand in 
u = (r ~ 2rhor)/''- This choice pushes the horizon to u ~ —1 and infinity to m = 1. The requirement that a series 
solution of the form ajU^ converges at both boundaries means that aj and '^j{—^yaj should both be finite; 

that is, that both a2j and a2j+i must be convergent. Substitution of the series ajU^ in the differential 
equation yields a five-term recurrence relation. This five-term recurrence relation can be split in two separate five-term 
recurrence relations, corresponding to the convergence conditions for even and odd coefficients. Finally, the even and 
odd five-term recurrence relations are reduced to three terms by two iterations of a Gaussian elimination step. 

A careful analysis of the extremal RN case is interesting for many reasons. First of all, as we have anticipated 
(see Figure |2l , the quasinormal mode spectrum for extremal RN black holes is characterized by an isospectrality 
between electromagnetic and gravitational perturbations: modes of with angular index I coalesce with modes of 
with index {I — 1) in the extremal limit. This has been shown to be a manifestation of supersymmetry 57.]. 
The supersymmetric nature of extremal Reissner-Nordstrom black holes is an attractive feature for the calculation of 
quantum effects, and it has been exploited in various contexts. Furthermore, topological arguments have been used to 
show that the entropy-area relation breaks down for extremal RN black holes jSg . So we believe that some caution is 
required in claiming that the connection between quasinormal modes and the area spectrum is still valid for extreme 
black holes, as recently advocated in [s^ . These problems may be connected with the finding that extremally charged 
black holes in a (non-asymptotically flat) anti-de Sitter spacetime could in principle be marginally unstable |60l [. 
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The explicit analysis of the extremal RN case confirms the above picture |4J|. It also shows that the standard 
Leaver method is surprisingly accurate even in this "extreme" situation, where it is not supposed to converge at all! 
Indeed, the largest error on the real frequency of the fundamental gravitational mode introduced by the "naive" Leaver 
treatment is only ^ 3 %. The authors of |4J| derived only a handful of quasinormal frequencies (namely the first 
three overtones for different values of I). As far as we know, results for extremal RN black holes in the intermediate 
damping regime have never been published to date. We show the resulting quasinormal spectrum in Figure (to be 
compared with Figure^. 

Extremal RN quasinormal frequencies have a similar pattern to Schwarzschild frequencies. Their real part first 
decreases, approaching the pure-imaginary axis as the overtone index grows. We find a quasinormal mode at uj = 
(0,-3.047876), while Chandrasekhar's formula 1)24(1 predicts a mode at w = (0,-3). Then ujji increases again, 
apparently approaching a constant value that could be compatible with the prediction (|27|l of the Motl-Neitzke 
asymptotic formula H25|) . Unfortunately it is hard to get stable numerical results for values of \ujj\ > 10 to explicitly 
verify this prediction. This relative numerical instability is not unexpected, since a) we need to solve simultaneously 
two different continued fraction relations, and b) for each of them we must carry out two Gaussian elimination steps. 



V. KERR BLACK HOLES 



A. Computational method 



As in the RN case, here we only give the basic steps of our computational procedure. More details can be found in 
[ioHi^ -ls H . For Kerr black holes, the perturbation problem can be reduced to two ordinary differential equations for 
the angular and radial parts of the perturbations, respectively. In Boyer-Lindquist coordinates, defining u = cos 9, 
the angular equation reads 



[(1 - u')Slm,u] „ + 



{auiu)'^ — 2aujsu + s +s Ai„i — 



(m + suY 



Sim — , 



(29) 



and the radial one 



ARlm,rr + (s + l)(2r - l)Rlm,r + V{r)Rim = , (30) 

where 

V{r) = {[{r'^ + a^fuj'^ ~ 2amujr + a'^m'^ + is{am{2r - 1) - uj{r'^ - a^))] 

+ [2isujr - a^u? A,™] } . (31) 

The parameter s = 0, —1, —2 for scalar, electromagnetic and gravitational perturbations respectively, a is the Kerr 
rotation parameter (0 < a < 1/2), and s-^im =s ^im(a'^) is an angular separation constant. In the Schwarzschild 
limit a = the angular separation constant can be determined analytically, and is given by the relation s^im = 
l{l + l)- sis + 1). 

Boundary conditions for the angular and radial equations yield two three-term continued fraction relations of the 
form Finding quasinormal frequencies is a two-step procedure: for assigned values of a, £, m and lu, first 

find the angular separation constant s^imiO''-^) looking for zeros of the angular continued fraction; then replace the 
corresponding eigenvalue into the radial continued fraction, and look for its zeros as a function of uj. In principle, 
the convergence of the procedure for modes with large imaginary parts can be improved, as described earlier, by a 
wise choice of the rest, Rn, oi the radial continued fraction. Expanding this rest as in formula (|11|) and introducing 
b = Vl — 4a2, we get for the first few coefficients: Co = —1, Ci = zLy/~2iujb, C2 = [3/4 -I- iuj{b + 1) — s]. 

The angular continued fraction, considered separately from the radial continued fraction, is a simple and efficient 
tool to compute the angular eigenvalues sAim{o-uj) of the spin- weighted spheroidal wave equations for complex values 
of the argument auj. We simply start from the known Schwarzschild limit, sAim{0) = l{l -I- 1) — s(s 4- 1). Then we 
solve the angular continued fraction, using this value as an initial guess, to compute s^im for generic values of acu. In 
section FV CI we show that this simple technique yields remarkably accurate results, even for large values of \auj\. 



B. Results 



The most relevant feature of the Kerr quasinormal spectrum is that rotation acts very much like an external 
magnetic field on the energy levels of an atom, causing a Zeeman splitting of the quasinormal mode spectrum. We 
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FIG. 7: "Zeeman-like" splitting of the fundamental gravitational mode with I = 2. We mark by dots the points corresponding 
toa = 0, 0.05,0.10, ...0.5. 





FIG. 8: Trajectory in the complex plane of the first eight Kerr quasinormal frequencies with m = 2 (solid lines) and m = —2 
(dashed lines). Filled circles mark the corresponding mode in the Schwarzschild limit. An arrow indicates the small loop 
described by the "exceptional" quasinormal mode with n = 6, that does not tend to the critical frequency for superradiance 
(see also Figures 3 and 4 in (gH ). 



illustrate the splitting of the fundamental quasinormal mode in Figure [7| (see also jSJ). As the rotation parameter a 
increases the branches with m = 2 and m = —2 move in opposite directions, being tangent to the the branches with 
m = 1 and m = — 1. In this simple case the effect of rotation on the modes is roughly proportional to m, as physical 
intuition would suggest. The mode with m = has a different behavior: it "moves" very quickly along a direction in 
the complex plane which is roughly perpendicular to the two branches with m ^ 0. 

In Figure |S1 we show the first eight quasinormal frequencies with m = 2 (solid lines) and m = —2 (dashed lines). 
For large enough damping, modes with m < are not tangent to the corresponding modes with positive m. Notice 
also that all modes with n < 9 and m > cluster at the critical frequency for superradiance ujsr = mfi, where fl is 
the angular velocity of the black hole horizon. This fact was first observed by Detweiler (bJ, who went on to suggest 
that all modes with positive m should tend to the critical frequency for superradiance as the black hole becomes 
extremal. It was even suggested that this could point to a marginal instability of extremal Kerr black holes. The 
study of the modes' excitation amplitudes in the same limit led to the conclusion that such an instability does not 
actually occur ^251, l63j . The mode with n = 6 (marked by an arrow) shows a peculiar deviation from the general 
trend, clearly illustrating the fact that some positive-m modes do not tend to the critical frequency for superradiance 
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FIG. 9: Left panel: real part of Kerr modes having m = as a function of a. Labels indicate the corresponding values of I and 
of the mode order n. Right panel; imaginary part of Kerr modes with / = 3, m = 0, n = 20 as a function of a. 



in the extremal limit. 

For gravitational wave emission, the dominant quasinormal mode is almost certainly the fundamental mode with 
I = m = 2. The reason is that a black hole formed in a stellar collapse or in a merger event is likely to have a rotating 
"bar" shape corresponding to spheroidal harmonic indices I = m = 2, and also that this mode is the most slowly 
damped, as we can see from Figure |H1 Echeverria [b^I found that the following analytic fit reproduces numerical data 
for the fundamental mode with I = m = 2 within about 10 %: 



Mujfi 



1-0.63(1 -a) 



3/10 



(0.37 + 0.19a), 



(32) 



~ (1.48 + 2.09a). (33) 

In these relations we restored the black hole masses; the last step is a series expansion for small values of the rotational 
parameter a. Fryer, Holz and Hughes produced a similar fit for the mode with I = 2, m = |65j| . 

In Figure |S1 we show the real and imaginary parts of some gravitational quasinormal frequencies with m = 0. The 
looping behavior is similar to that of RN modes modes, and the number of loops increases with the damping of the 
mode. This same looping behavior was recently found by Glampedakis and Andersson for scalar perturbations of 
Kerr black holes, using a different method [ssj . 

Figure [TUI shows the trajectory in the complex plane of some scalar modes with m = 0. The spirals wind up very 
quickly as the damping grows, and the "center" of the spiral rapidly moves towards the pure-imaginary axis. In 
section FVPI we will see that the oscillation frequency of modes with m = probably tends to zero as \ojj\ —> oo. 

In [^i we found that for intermediate damping (say 10 < n < 50) quasinormal frequencies belong to three different 
families depending on whether m > 0, m = or m < 0. Figure lTTl shows some representative results for electromagnetic 
modes (s = ^l)i but the qualitative behavior of gravitational and scalar modes is the same: 

1) Modes with m > 0. For 10 < < 50, tj/j has a relative minimum as a function of a, but the numerics do not 
show any convergence to an asymptotic limit. Only the real part of gravitational modes with I = m = 2 tends 
to the limit ujr = 2Q when n ~ 50 (see and Figure 1 in This deceiving convergence seems related to 
the peculiar behavior of the angular separation constant for s\ — I — m — 2 (see Figure [T^ . For m > most 
modes in the range 10 ^ n < 50 approach the limiting value ajR = mfl as a —^ 1/2, but there are exceptions (see 
Figure IHl). The imaginary parts of all modes with m > have a separation 2ttTh, where now Th is the Hawking 
temperature of a Kerr black hole. The agreement with the separation of consecutive Schwarzschild quasinormal 
modes does not seem to be a coincidence, Tu being a non-trivial function of the rotation rate. 

2) Modes with m = 0. The real and imaginary parts of these modes oscillate as a function of a, a behavior 
reminiscent of RN quasinormal modes. As n increases, (jJr becomes smaller. The spacing between the imaginary 
parts oscillates and is not given by 2'kTh- 
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FIG. 10: Trajectories of a few scalar modes with I — m — 0. The different panels correspond to the fundamental mode (top 
left), which does not show a spiralling behavior, and to modes with overtone indices n — 2, 4, 10. 



3) Modes with m < 0. At least for a > 0.25, loj^ seems to converge to a (weakly a-dependent) limit lor ~ mzu, 
where -cu ~ 0.12, whatever the value of I and the spin of the perturbing field. The spacing in the imaginary 
parts for these values of n does not show convergence, and is not given by 2t:Th- 

These results are even more puzzling than in the RN case. To probe the asymptotic regime we clearly need to push 
the calculation to higher damping. The main problem here is that Leaver's approach requires the simultaneous solution 
of the radial and angular continued fraction conditions. For mode order n > 50 the method becomes unreliable, or it 
just fails to provide results. A way around this "coupling problem" is to study the asymptotic behavior of the angular 
equation as \aLu\ ~ \aLUi\ — > oo. In the next section we carry out this task. We first compute numerically s-^imio-^) 
when uj ~ iw/ and |acij| —>■ oo. Then we determine analytically the leading-order behavior of the separation constant. 
Finally, in section lV Dl we replace this analytic expansion in the radial continued fraction, effectively getting rid of the 
angular equation. By this trick we can find quasinormal frequencies with n ^ 50 solving only the radial continued 
fraction, and ultimately determine the asymptotic behavior of ujji{a) for Kerr black holes. 



C. Angular eigenvalues 

The anal ytic al p roperties of the angular equation (|29|1 and of its eigenvalues have been studied by many authors 
lil; |6E 120; llE lzl> Izll ' Series expansions of s^Zm for \auj\ <C 1 have long been available, and they agree well 
with numerical results (see eg. [tJ, where some mistakes and notational differences in the previous literature are 
discussed and corrected). 
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Series expansions for large |aa;| have been known for a fong time in the special case of s = 0. Such an analytical 
power-series expansion for large (pure real and pure imaginary) values of ato can be found, for example, in Flammer's 
book fpiil- Flammer's results are in good agreement with the pioneering numerical work by Oguchi '73|, who computed 
angular eigenvalues for complex values of aw and s = 0. A review of numerical methods to compute eigenvalues and 
eigenfunctions for s — can be found in l67|. Recentlv Barrowes et al. studied the large frequency behavior of s = 
harmonics in the complex frequency plane |68|. while Falloon et al. developed a Mathematica notebook to compute 
s = harmonics for general complex values of the frequency . 

Unfortunately similar results for large values of |aw| and general spin s are lacking. Even worse, until recently the 
few analytical predictions were contradictory M- The situation fo^ l^rgereal frequencies has finally been clarified 
by Casals and Ottewill [t^, who corrected some mistakes in reference [T^- Here we fill the remaining gap. We 
present numerical and analytical results for s^im when the frequency acu is large and purely imaginary. In Flammer's 
language we are looking at the prolate case, while Casals and Ottewill deal with the oblate case. 

We carry out our numerical calculations choosing acu to be pure imaginary and solving Leaver's angular continued 
fraction as explained in section FV Al Our numerical values for s — and purely imaginary frequency (top left panel 
in Figure IT^ are in agreement with Tables 10-12 in 66] . For complex oj and s = we could reproduce Table 2 in 
[7J|. As a final sanity check, we verified in a few representative cases that the Mathematica notebooks presented in 
|67l l69| are in agreement with our code when we specialize to s = 0. 

We also verified numerically some elementary analytical properties of the eigenvalues: for example, the computed 
eigenvalues satisfy the relation _s^!m =s + 2s. In the oblate case considered in [t^I the eigenvalues grow 
quadratically, s^Zm ~ "l^'^Pj as they shoul d 1761 . Preliminary comparisons with Casals show excellent agreement on 
the computed eigenvalues for all values of s |77|. 
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FIG. 12: Top left: angular separation constant aAim for s = as a function of c = au>i. The numerical behavior is in 
agreement with Flammer's asymptotic formula for the prolate case, and the separation constant is purely real. The slopes kim 
c are either 1 or 3. Top right, bottom left: angular separation constant for s — ~1 and s — —2. Now a Aim is 
not purely real, but because of eq. 13511 it is sufficient to compute eigenvalues with m > 0. In both plots, lines starting from 
aAim = l{l + 1) — s(s + 1) are the real parts, and lines starting from aAim = are the imaginary parts. Notice the "kink" in 
the real part at c ~ 3 when s = —2 and I = m = 2: there is no discontinuity, just a change of slope. This peculiar behavior 
might explain the deceiving convergence of the quasinormal frequency to lor = 2f2 that we observed in Q for I — m — 2 and 
s — —2. Bottom right: real parts of aAim for s = —1 (in black) and s = —2 (in red). The slopes for the linear growth at large 
c are all consistent with Equation I38II . 



For s = the series expansion given by formula (8.1.11) in |6a | is an excellent approximation to the numerical data: 

qAi^ = (2L + l)c-(2L2 + 2i + 3-4m2)2-2-(2L + l)(i2 + i + 3_8TO2)2-''c-^ (34) 
- [5(L'' + + 7L + 3) - \%m^{2l? + 2L+ l)]2-^c-'^ + + ©(c^^) . 

Here we defined c — aui and L = I — \m\. Comparison with numerical results (see the top left panel of Figure [T^ 
shows that this expansion is reasonably accurate down to c ~ 2. 

For s — ~1 and s = —2 the angular separation constant is complex. We can limit our calculations to positive m's 
because of the following symmetry property (see eg. |40|): 

sj4; _„j =s j4*„ , (35) 

where the asterisk denotes complex conjugation. Numerical results for the lowest radiatable multipoles are given in 
Figure 1121 There we plot together real and imaginary parts of the eigenvalues. It is apparent that the real part is 
dominant, growing linearly as c ^ oo. 
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The proportionality constant for this hnear growth can be found by a straightforward generahzation of Flammer's 
method to general values of s 0. Define a new angular wavefunction Zim{u) through |7ll | 

= (l-u2)^Zi™(M), (36) 
and change independent variable by defining x — ^/2cu^ where c? = —[au)'^. Substitute this in 129|) to get: 



4, 



iylcx — m{m + 1) 



2msx 



Zim + {2c~X^)Zir, 



2{m + s + l)xZim,x = . 



(37) 



2 ^ ' \/2c + xi 

When c ^ oo, this equation becomes a parabolic cylinder function. The arguments presented in [s^ l7ll lead to 

= (2L + l)c + 0(c") , c->oo, (38) 
where L is the number of zeros of the angular wavefunction inside the domain. One can show that 



L 



I - \m\ , |m| > |s|, 
l-\s\, \m\<\s\. 



(39) 



Higher order corrections in the asymptotic expansion can be obtained as indicated in pq . We will not need them in 
our calculation of highly damped modes. 



D. High damping 

To compute the asymptotic quasinormal frequencies of the Kerr black hole we use a technique similar to that 
described in |4^ . We fix a value of the rotation parameter a. We first compute quasinormal frequencies for which 
I aw I 1, so that formula H38I) is only marginally valid. This procedure is consistent with our previous intermediate- 
damping calculations: for example, when we include terms up to order |act;|^^ in the asymptotic expansion for oAim 
provided in [g^ . our new results for a ~ 0.1 and I — m = 2 match the results for the scalar case at overtone numbers 
20 < n < 30. Then we increase the overtone index n (progressively increasing the inversion index of our "decoupled" 
continued fraction). Finally, we fit our numerical results by functional relations of the form: 

ujR{n) LJR + + J^^a^ + Lu'^R^a^, (40) 

where a — 1/|ij-'/| or a — -y/T/jwTf. At variance with the non-rotating case fits in powers of l/jw/l perform better, 
especially for small and large a. However, both fits break down as a ^ 0: the values of higher-order fitting coefficients 
increase in this limit, so that subdominant terms become as important as the leading order, and the extraction of 
the asymptotic frequency lur becomes problematic. The numerical behavior of subdominant coefficients supports the 
expectation (which has not yet been verified analytically) that subdominant corrections are a-dependent. Therefore 
one has to be careful to the order in which the limits n ^ oo, a are taken 0, 0, These observations are 
consistent with the fact that, in the RN case, the zero-charge limit of the asymptotic quasinormal frequency spectrum 
does not yield the asymptotic Schwarzschild quasinormal spectrum. 

We have extracted asymptotic frequencies using two independent numerical codes. For each value of a, we found 
that the extrapolated value of ujr is independent of s (Fig. I13|l . independent of I and proportional to m (Fig. I14|l : 

lur, ~ mzu{a) . (41) 

We obtained lur computing quasinormal frequencies both for scalar perturbations (s = 0) and for gravitational 
perturbations (s = —2). In both cases we picked I ~ m — 2. The agreement between the extrapolated behaviors of 
ujR as a function of a is excellent, suggesting that both sets of results are typically reliable with an error < 1%. Our 
results are also weakly dependent on the number of terms used in the asymptotic expansion of s-Aim- this provides 
another powerful consistency check. We have tried to fit the resulting "universal function" , displayed in Fig. ^1 by 
simple polynomials in the BH's Hawking temperature T and angular velocity 11 (and their inverses). None of these 
fits reproduces our numbers with satisfactory accuracy. It is quite likely that asymptotic quasinormal frequencies 
will be given by an implicit formula involving the exponential of the Kerr black hole temperature, as in the RN case 

For any a, the imaginary part ui grows without bound. Quite surprisingly, the spacing between modes Suji is a 
monotonicall y in creasing function of a: it is not simply given by 2ttTh, as some numerical and analytical calculations 
suggest A power fit in a of our numerical results yields: 

Slui = 1/2 + 0.0438a - 0.0356a^ (42) 
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FIG. 13: Real part of high-order quasinormal frequencies for scalar (s = 0), electromagnetic {s = —1) and gravitational 
(s = —2) perturbations of a Kerr BH with a — 0.1 (I = m = 2). Quasinormal frequencies of different spins converge to the 
same value. For any kind of perturbation we are already deep in the region of validity of the asymptotic expansion (1381 . since 
I act; I ~ 60. 




FIG. 14: Gravitational quasinormal frequencies of a Kerr BH with a = 0.1, I = 2 and different values of m. Dashed lines and 
open circles correspond to negative m. The asymptotic value does not depend on the sign of m, and is proportional to \m\. 
On the right we show m = Q quasinormal frequencies for different values of the rotation parameter: the larger a, the faster ujr, 
tends to zero as |a;j| oo. 



VI. ALGEBRAICALLY SPECIAL MODES 



Algebraically special modes of Schwarzschild black holes have been studied for a long time, but only recently their 
understanding has reached a satisfactory level. Among the early studies rank those of Wald [Tg and of Chandrasekhar 
ji^ . who gave the exact solution of the Regge- Wheeler, Zerilli and Teukolsky equations at the algebraically special 
frequency. The nature of the quasinormal mode boundary conditions at the Schwarzschild algebraically special 
frequency is extremely subtle, and has been studied in detail by Maassen van den Brink [79i |. 

To understand the problem, it is useful to recall that black hole oscillation modes can be classified into three groups 
(here we use Maassen van den Brink's "observer-centered definition" of the boundary conditions): 

1) "standard" quasinormal modes, which have outgoing wave boundary conditions at both sides (that is, they are 
outgoing at infinity and "outgoing into the horizon" ) ; 

2) total transmission modes from the left (TTM^'s) are incoming from the left (the black hole horizon) and outgoing 
to the other side (spatial infinity); 
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FIG. 15; Asymptotic real part ur, — 2zu{a) of the I = m — 2 gravitational and scalar quasinormal frequencies extrapolated 
from numerical data: ljr 2cc7(l/2) ~ 0.21 as a 1/2. Results are independent of I, s and the sign of m. 



3) total transmission modes from the right (TTM^j's) are incoming from the right and outgoing to the other side. 

The Schwarzschild "algebraically special" frequency is given by formula (|13|l for Schwarzschild black holes [or by 
formula H24|) for RN black holes], and traditionally it has been associated with TTM's. In fact, the algebraically 
special frequency is so called precisely because of the requirement that the perturbations should be "special": they 
should contain only ingoing waves (the perturbed Weyl scalar 'i'o ^ 0) or only outgoing waves (^'4 ^ 0). Speciality is 
equivalent to the condition that a certain quantity (the Starobinsky constant) should be zero. When Chandrasekhar 
determined the perturbation satisfying this condition pal he did not check that it does indeed satisfy TTM boundary 
conditions. It turns out that, in general, it does not [73 ■ The Regge- Wheeler equation and the Zerilli equation (which 
are known to yield the same quasinormal mode spectrum, being related by a supersymmetry transformation) have to 
be treated on different footing at fli , since the supersymmetry transformation leading to the proof of isospectrality 
is singular there. In particular, the Regge- Wheeler equation has no modes at all at fli, while the Zerilli equation has 
both a quasinormal mode and a TTM^. 

Numerical calculations of algebraically special modes have yielded some puzzling results. Studying the Regge- 
Wheeler equation (that should have no quasinormal modes at all according to Maassen van den Brink's analysis) and 
not the Zerilli equation. Leaver 0| found a quasinormal mode which is very close, but not exactly located at, the 
algebraically special frequency (Table QJ. Namely, he found quasinormal modes at frequencies fij such that 

Cl'^ = 0.000000 - 3.9980001, fig = -0.000259 - 20.0156531. (43) 

Notice that the "special" quasinormal modes fij are such that 3? iil2 < \^2\, ^ 1^3 > l^sl, and that the real part of 
rig is not zero. Similarly, in the extremal RN case we find a quasinormal frequency that is very close, but not exactly 
equal to, the algebraically special frequency found by Chandrasekhar. Maassen van den Brink 't^ speculated that 
the numerical calculations may be inaccurate and that no conclusion can be drawn on the coincidence of fli and fij, 
"if the latter does exist at all". 

An independent calculation was carried out by Andersson Is^. Using a phase-integral method, he found that the 
Regge- Wheeler equation has pure imaginary TTMi^'s which are very close to fti for 2 < Z < 6. He therefore suggested 
that the modes he found coincide with tli, which would then be a TTM. Maassen van den Brink observed again 
that, if all figures in the computed modes are significant, the coincidence of TTM's and quasinormal modes is not 
confirmed by this calculation, since ilj and £7; are numerically (slightly) different. 

Onozawa showed that the Kerr mode with overtone index n — 9 tends to il; as a ^ 0, but suggested that 
modes approaching (li from the left and the right may cancel each other at a = 0, leaving only the special (TTM) 
mode. He also calculated this (TTM) special mode for Kerr black holes, solving the relevant condition that the 
Starobinsky constant should be zero and findin g th e angular separation constant by a continued fraction method; his 
results improved upon the accuracy of those in '46'. 

The analytical approach adopted in j79l | clarified many aspects of the problem for Schwarzschild black holes, but the 
situation concerning Kerr modes branching from the algebraically special Schwarzchild mode is still far from clear. In 
|79|| Maassen van den Brink, using slow-rotation expansions of the perturbation equations, drew two basic conclusions 
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FIG. 16: The left panel shows the trajectories described in the complex-tj plane by the doublets emerging close to the 
Schwarzschild algebraically special frequency {CI2 ~ —4i) when m > and 1 — 2. Notice that the real part of modes with 
m > tends to uir = m as a —> 1/2. The right panel shows the spiralling trajectory of the mode with m = 0. 



on these modes. The first is that, for a > 0, the so-called Kerr special modes (that is, solutions to the condition 
that the Starobinsky constant should be zero [46ll6ll |) are all TTM's (left or right, depending on the sign of s). The 
TTM/j's cannot survive as a — > 0, since they do not exist in the Schwarzschild limit; this is related to the limit a — > 
being a very tricky one. In particular, in this limit, the special Kerr mode becomes a TTM^ for s = —2; furthermore, 
the special mode and the TTMi,; cancel each other for s = 2. Studying the limit a ^ in detail, Maassen van den 
Brink reached a second conclusion: the Schwarzschild special frequency f2; should be a limit point for a multiplet of 
"standard" Kerr quasinormal modes, which for small a are well approximated by 

^. 33078176 3492608. , 4, 

when / = 2, and by a more complicated formula - his equation (7.33) - when I > 2. We numerically found quasinormal 
modes close to the algebraically special frequency, but they do not agree with this analytic prediction (not even when 
the rotation rate a is small). 

Maassen van den Brink suggested (see note [46] in 't^) that quasinormal modes corresponding to the algebraically 
special frequency with m > may have one of the following three behaviors in the Schwarzschild limit: they may 
merge with those having to < at a frequency f2J such that \Ct[\ < (but \Ct'i\ > for I > 4) and disappear, 
as suggested by Onozawa |6lj : they may terminate at some (finite) small a; or, finally, they may disappear towards 
u) — —100. Recently Maassen van den Brink et al. (sj put forward another alternative: studying the branch cut on 
the imaginary axis, they found that in the Schwarzschild case a pair of "unconventional damped modes" should exist. 
For I = 2 these modes were identified by a fitting procedure to be located on the unphysical sheet lying behind the 
branch cut (hence the name "unconventional") at 

w± = T0.027+ (0.0033- 4)1. (45) 

An approximate analytical calculation confirmed the presence of these modes, yielding 

UJ+ ~ -0.03248 + (0.003436 - 4)i, (46) 

in reasonable agreement with (|45|l . If their prediction is true, an additional quasinormal mode multiplet should 
emerge near J7; as a increases. This multiplet "may well be due to uj± splitting (since spherical symmetry is broken) 
and moving through the negative imaginary axis as a is tuned" In the following paragraph we will show that a 
careful numerical search indeed reveals the emergence of such multiplets, but these do not seem to behave exactly as 
predicted in [sj]. 
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FIG. 17: The top row shows the real and imaginary parts (left and right, respectively) of the "doublet" of quasinormal modes 
emerging from the algebraically special frequency as functions of a. The doublets only appear when m > 0. We also overplot 
the real and imaginary parts of the mode with Z = 2, m. = (showing the usual oscillatory behavior). The bottom row shows, 
for completeness, the real and imaginary parts (left and right, respectively) of modes with negative m branching from the 
algebraically special frequency. 



A. Numerical search and quasinormal mode multiplets 



As we have seen from the summary in the previous paragraph, the situation for Kerr modes branching from the 
algebraically special Schwarzschild mode is not clear. There are still many open questions. Is a multiplct of modes 
emerging from the algebraically special frequency when a > 0? Can quasinormal modes be matched by the analytical 
prediction (|44|l at small values of a? If a doublet does indeed appear, as recently suggested in .81], does it tend to 
the algebraically special frequency ^i = — 4i as a — > 0, does it tend to the values predicted by formula H45|l . or does 
it tend to some other limit? 

After carrying out an extensive numerical search, we have found some surprises. Our main result is shown in the 
left panel of Figure [TCI There we show the trajectories in the complex plane of quasinormal modes with / = 2 and 
m > 0: a doublet of modes does indeed appear close to the algebraically special frequency. Both modes in the doublet 
tend to the usual limit {(I2 = m) as a ^ 1/2. We have tried to match these "twin" modes with the predictions of 
the analytical formula H44|l . Unfortunately, none of the two branches we find seems to agree with l)44|l at small a. We 
could only find a mode doublet for m > 0. For m < the behavior of the modes is, in a way, more conventional. 
For example, in the right panel of Figure^jwe see the I — 2, m = mode emerging from the standard algebraically 
special frequency CI2 and finally describing the "usual" spirals as a increases. 

In the top left panel of Figure El we see that the real part of the twin modes with m > goes to zero as a ^ with 
an m-dependent slope. However, the top right panel in the same Figure shows that the imaginary part of the modes 
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does not tend to —4 as a ^ 0. Qualitatively this behavior agrees rather well with that predicted by equation H45|l . 
Extrapolating our numerical data to the limit a — *■ yields, however, slightly different numbers. Our extrapolated 
values for ? = 2 are u = (-4 - 0.10)i and uj = (-4 + 0.09)i. 

At present, we don't know why the doublet only appears when to > 0. This fact is confirmed by numerical searches 
close to the algebraically special frequency i^a for / — 3. Once again, a quasinormal mode multiplet only appears when 
TO > 0. In particular, we see the appearance of a doublet similar to the modes shown in the left panel of Figure ITHl 
Extrapolating numerical data for the I = 3 doublet in the limit a ^ yields w = (—20 — 0.19)i and w = (—20 + 0.24)1. 

A more careful search near the algebraically special frequency fl^ surprisingly revealed the existence of other quasi- 
normal modes. However, the additional modes we find may well be "spurious" modes due to numerical inaccuracies, 
since we are using the continued fraction technique in a regime S> 1, wn ~ 0) where it is not supposed to work 
well. 



VII. CONCLUSIONS AND OUTLOOK 

The quasinormal mode spectrum of Schwarzschild, RN and Kerr black holes is overwhelmingly rich and complex. 
Our understanding has certainly improved over time, but many mysteries are still unsolved. The following is a (partial 
and necessarily biased) list of open problems. 

1) A general lesson from the study of charged and rotating black holes seems to be: the spectrum of the limit is 
not the limit of the spectrum. This is a well-known fact in spectral theory, but it is quite striking in the black 
hole context. "Ordinary" black hole quasinormal modes suggest that there should be some continuity between 
the Schwarzschild and Kerr (RN) solutions as a ^ (Q ^ 0). At a closer look, the asymptotic spectrum in the 
limit — > oo seems to violate this continuity requirement. What is the mathematical and (more importantly) 
physical motivation of this discontinuity? 

2) How does the picture change when we include both rotation and charge? As anticipated in the introduction, the 
main difhculty here is that perturbation equations for Kerr-Newman black holes are non-separable (unless we 
introduce some ad hoc approximations). Work in this direction is ongoing '84'|. 

3) The algebraically special mode is a mystery on its own. Numerical methods show the presence of a quasinormal 
mode close to (but not quite at) the algebraically special frequencies determined by Chandrasekhar ~ this is 
true both in the Schwarzschild and in the extremal RN cases (see Figures ^ and O . In the Kerr case, a doublet 
of modes with to > comes out of the algebraically special Schwarzschild frequency, but no mode splitting is 
observed when m < 0. Analytical work on the "supersymmetry breaking" occurring at the algebraically special 
frequency only partially clarifies the situation [79l |. We definitely need a better understanding of the meaning 
of quasinormal modes and total transmission modes in this situation. 

4) The behavior of Kerr quasinormal frequencies with m > is now better understood. Not all modes tend to the 
critical frequency for superradiance in the limit a ^ 1/2. Indeed, for some modes cuj does not even tend to zero 
in the extremal limit (see Figures |S1 and I15() . The explanation seems to be that some implicit assumptions in 
Detweiler's original treatment were probably overlooked ,62^ . 

5) Our present physical understanding of highly-damped black hole oscillations is admittedly very poor. What 
is the meaning of the spirals we see for RN quasinormal modes (and for Kerr modes with m = 0)? Is the 
"universal" Kerr asymptotic frequency zu(a) (Figure^)) somehow connected to other fundamental frequencies 
of the Kerr spacetime, eg. last stable photon orbits? Why does the RN asymptotic spectrum depend on the 
causally disconnected region inside the horizon? Is the same true of the Kerr spacetime? And what about 
Kerr-Newman black holes? 

At a more fundamental level, we should go back to the original motivations behind this investigation. The complexity 
of the (purely classical) black hole oscillation spectrum is striking. Can we still believe that there is a link between 
asymptotic quasinormal frequencies and black hole quantization, or was it just a numerical coincidence? Can we 
really think of black holes as "quantum gravity atoms" 0,^3, being the simplest objects we can build out of gravity 
alone? If the analogy between classical black hole oscillations and atomic spectra holds, the bizarre nature of the 
quasinormal mode spectrum is telling us that there's a long way to go before we fully understand how to quantize 
gravity. 
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